XY models with disorder and symmetry-breaking fields in two dimensions 
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The combined effect of disorder and symmetry-breaking fields on the two-dimensional XY model 
is examined. The study includes disorder in the interaction among spins in the form of random phase 
shifts as well as disorder in the local orientation of the field. The phase diagrams are determined and 
the properties of the various phases and phase transitions are calculated. We use a renormalization 
group approach in the Coulomb gas representation of the model. Our results differ from those 
obtained for special cases in previous works. In particular, we find a changed topology of the phase 
diagram that is composed of phases with long-range order, quasi-long-range order, and short-range 
order. The discrepancies can be ascribed to a breakdown of the fugacity expansion in the Coulomb 
gas representation. Implications for physical systems such as planar Josephson junctions and the 
faceting of crystal surfaces are discussed. 

PACS numbers: 64.60.Ak; 75.10.Nr 
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I. INTRODUCTION 

We reconsider the two-dimensional XY model with 
random phase shifts and a symmetry-breaking field, 
which is described by the reduced Hamiltonian 

U = -K a ^cos(V a 6 r ~ A r>a ) - H o ^2cos(p0 r - 4> T ). 

r,a r 

(1) 

XY spins 9 r are placed on the sites r of a square lat- 
tice. The coupling of nearest neighbors with a spin stiff- 
ness K = J/T favors a relative angle prescribed by the 
quenched random phase shifts A, which we assume to 
be Gaussian distributed with zero mean and variance 
A 2 a = a. The global rotation symmetry of spins is bro- 
ken by the field Hq that selects p £ N equivalent favorable 
directions for the spins. The model will be considered 
with uniform fields (4> = 0) and with random local orien- 
tations 4> °f the field. 

This model describes a variety of physical systems, 
ranging from disordered magnetstJ to Josephson-junction 
arraysa and planar Josephson-junctions.B It is also closely 
related to surfaces of crystals with bulk disorderu and 
even to periodic elastic media such as-,vortex lattices in 
the presence of a pinning potential,BEl two-dimensional 
crystals with quenched impurities!] or adsorbates on dis- 
ordered substratesBB 

The study of the above and related models started al- 
ready decades ago. We briefly recall some marked steps 
in its history to provide the background for the present 
work. The XY model in the absence of disorder and mag- 
netic fields exhibits the Kosterlitz-Thouless (KT) tran- 
sition from a phase with quasi-long-range order (QLRO, 
i.e., algebraic decay of the spin-spin correlation) to short- 
range order (SRO, i.e., expcmeptial decay of the spin-spin 
correlation) at = 2/7r,t2r£3 where the effective spin 



stiffness Koo < Kq is renormalized by thermal fluctua- 
tions. The transition was identified with an unbinding of 
vortex pairs. Weak uniform fields have been-included by 
Jose, Kadanoff, Kirkpatrick, and Nelson.E3Ej For p > 4 
they found that the phase with QLRO is stable against 
such fields over a temperature range 8tt/p 2 < K^ 1 < n/2 
and that weak fields become relevant at low tempera- 
tures < 8w/p 2 where they stabilize long-range order 
(LRO, i.e. finite magnetization). For p < 4 the transi- 
tion resulting from the competition between the ordering 
tendency of the magnetic field and the disordering ten- 
dency of the thermal fluctuations could not be captured 
by a KT-like approach. 

Houghton, Kenway, and YingtS found the phase with 
QLRO to be unstable to magnetic fields with random 
orientation only for K^ 1 < An/p 2 . Accordingly, a phase 
with QLRO can be stable within a finite temperature 
range in/p 2 < K^ 1 < it/ 2 provided p > 2^/2- Their 
study was extended by Cardy and OstlundE£l (CO) to 
address the nature of the phase transitions and of the 
low-tempcrature phase. Their analysis implied that ran- 
domness of the direction of the magnetic field effectively 
generates a new type of disorder on large scales (in the 
sense of jsurenormalization group), namely random phase 
shifts AJB 

Rubinstein, Shraiman, and NelsonEJ (RSN) examined 
the XY model with random phase shifts in the absence 
of fields. They found QLRO to be stable against dis- 
order below a certain critical strength and a phase dia- 
gram with reentrant shape. Paczuski and KardarEJ (PK) 
added uniform fields and obtained a phase diagram com- 
posed of various phases with LRO, QLRO, and SRO. 
The transitions between these phases showed multiple 
reentrance. 

Recent advance in the understanding of the bond- 
disordered model without fields has led to a revised non- 
reentrant phase diagram, where the phase with QLRO 
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turned out to be more stable than found by RSN.I 
The earlier underestimation of the order in the system 
was caused by a subtle breakdown of the fugacity expan- 



sion in the Coulo: 
low temperatures ' 24 



as representation of the model at 
This breakdown was overcome in 
these recent approaches by nonperturbative techniques. 
The features of the revised phase jdiagram have been con- 
firmed by numerical simulations. Ej'EZI 

The breakdown of the fugacity expansion is likely to 
affect all problems involving random phase shifts. In par- 
ticular one has to ask how the structure of the phase di- 
agram, the nature of the phases and the nature of the 
phase transitions are modified. It is also important to 
examine whether a similar breakdown occurs in a per- 
turbative treatment of the fields. The purpose of the 
present work is to reexamine model (Q) in view of these 
questions. 

The method of our choice is a mapping of model (|l|) 
on two coupled Coulomb gases, "vortices" and "charges" 
(Sec. II). Thereby we take a pre-averaging over disor- 
der using the replica trick. We establish renormalization 
group flow equations for the model parameters to cap- 
ture the large scale properties of the system (Sec. III). 
From the evaluation of these flow equations (Sec. IV) we 
obtain results that are significantly different from those 
obtained by CO and PK for the special cases of ferro- 
magnetic bonds with random fields and bond disorder 
with uniform fields. Our phase diagrams are illustrated 
in Fig. ^ for these two previously examined special cases. 
We conclude with a discussion of our results and of con- 
sequences for physical realizations of the XY model (Sec. 
V). Appendices with some calculational details and a dic- 
tionary for a translation of notations from previous works 
to the present work are included. 



II. MODEL TRANSFORMATIONS 

In this Section we specify the stochastic properties of 
the random phase shifts considered subsequently. In ad- 
dition, we transfer randomness in the orientation of the 
symmetry-breaking fields to the random phase shifts such 
that we need to deal only with (correlated) bond disor- 
der. The latter model is then mapped onto two coupled 
Coulomb gases. Eventually we arrive at a replica repre- 
sentation of the model that serves as basis for our further 
analysis. 



A. Disorder specification 

In order to be able to treat the model irrespectively of 
a randomness in the local field orientation, we decompose 
the field A = {A a } in Fourier space into its longitudinal 
and transverse component, 



(2b) 



Al, a =^P k ^, /3 A ki , = £ ( i _ ] lk: , (2( , 
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Both components (polarizations II = L, T) are supposed 
to be Gaussian distributed with zero average and to have 
variances (over-lining denotes averaging over quenched 
disorder) 



^;A=^k%<5(k + k'), 



(3) 



where we explicitly allow for <r L ^ <r T . Then the random 
phase shifts have a variance a = (a L + tr T )/2 and are 
in general correlated in space. Only in the special case 
a L = a T the phase shifts are uncorrelated on different 
bonds of the lattice. 

Neglecting the effect of vortices and fields completely, 
the spin-spin correlation (angular brackets denote a ther- 
mal average for a fixed disorder realization) 



T p (r) = (cos\p8 r - p6 ]) 



(4) 



can be calculated in the spin-wave approximation and is 
found to decay algebraically, 



r P (r) 



(5) 



with an exponent r\ = (p 2 /2ir)(l/K + a L ). In this ap- 
proximation 77 depends only on the longitudinal disorder. 
As in Ref. we will also examine correlations of the 
Edwards- Anderson order parameter 



e 



iei-iei 



{ij; ab (r)V> ab (0)) - r"" 



(6a) 
(6b) 



between different replicas a =/= b of the system. In the 
spin- wave approximation the exponent rj = 1 /irK is com- 
pletely independent of disorder £3 



B. Random fields 

A convenient way of dealing with the randomness in 
the local field orientation </> is to eliminate the angles <f> 
in favor of modified phase shifts A = A V0/p by a 
substitution 9 = 9 + <f>/p . This transformation leaves 
the thermodynamic partition sum invariant, i.e., 



Z({A},{0})= / VOe 



-H 



Z({A},{0}). 



(7) 



A = A 1 



(2a) 



Therefore, the existence and location of phase transitions 
in the phase diagram are not affected by this manipula- 
tion. In contrast, spin correlation functions are not in- 
variant and can display qualitatively different large-scale 
behavior before and after the transformation. 
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In order to restrict our analysis of random phase shifts 
to those with correlations (0) even after the transforma- 
tion, we limit the possible distributions of the angles <fi 
to those having a Gaussian distribution with correlations 



Mv = A— S(k + k'). 



(8) 



Only such random orientations are transformed by the 
above mapping into longitudinal phase shifts with a vari- 
ance 



(9) 



Equation (g) implies correlations (4> r — <fio) 2 ~ (A/V) lnr 
for r > 1 in real space. However, in the limit A — > oo 
the relative fluctuations even of neighboring angles are 
that large that they can be considered practically as un- 
corrected since they enter the Hamiltonian only modulo 
2n. In this sense the examination of the model ([j]) with 
random phase shifts and uniform field effectively includes 
the model with random symmetry-breaking fields in the 
limit a L — > oo. 



C. Coulomb gas representation 

We apply standard techniquesEl to map model ([!]) onto 
a Coulomb gas_.For this purpose we first use the Villain 
approximation!^ to replace both cosines in Eq. (Q) by a 
slightly different periodic potential and obtain an effec- 
tive Hamiltonian 



w = f £(v< 

r 

H 



27rM r 



— ^(p6» r + 27rm r ) 2 



(10) 



The partition sum of this effective model contains not 
only the integration over the angles 9 but also a sum- 
mation over the integer-valued vector field M and the 
integer-valued scalar field m. The parameters of this ef- 
fective model are related to the original parameters .by 
K = f(K ) and H = f(H ) through the functionHB " 



fix) 



T/lm> 2 /4) 



for 
for 



x > 1, 
x < 1. 



(11) 



Since we are mainly interested in the stability of the 
phase with QLRO (i.e., Kq > 1) to symmetry-breaking 
fields (i.e., Hq <C 1), we effectively operate in the limiting 
cases K « K and H w -l/ln(i$/4). 

For small H the sum over m converges very slowly. 
Therefore it is of advantage to go over from the field m to 
a conjugate integer- valued field n by means of the Poisson 
summation formula (some calculational details are given 
in Appendix |A|) . After integrating out the angles 9 in the 
partition sum we end up with a new effective Hamiltonian 



for two coupled Coulomb gases. The first Coulomb gas 
is the gas of vortices with a vorticity JV = V x M 6 Z. 
These vortices are located on the dual lattice, i.e., on 
plaquette centers The "charges" n r eZ constitute 
the second Coulomb gas. These Coulomb gases are cou- 
pled to "quenched vortices" Qr = (1/2tt)V x A and 
"quenched charges" q r = (i^/p)V • A, which are gener- 
ated by the transverse and longitudinal components of 
the bond disorder respectively. 

The Hamiltonian of the effective model reads 



tL — Ti. v + Ti c + Tlvc + n-x-, 



(12a) 



U, = - y - V ln|R-R'|(JVR-QR)(JVR,-Q R 

R^R' r 

+E V J2K, (12b) 

R 

Ho = -- ^^ ln l r - r 'IK -%)K' -iq T >) 



r 

n vc =zp^w(R-r)A R 7 



R,r 



<7r<7r 



(12c) 
(12d) 
(12e) 



It has been split into the interaction among vortices Tt v , 
the interaction among charges 7i c , the interaction be- 
tween vortices and charges H vc , and an excess term H x . 
Particles of the same type interact logarithmically and 
the coupling between the two gases is given by the angle 



w(R — r) = arctan(R — r) 



(13) 



enclosed by their relative distance and an arbitrary ref- 
erence direction in the xy plane (see Appendix H for 
some intermediate steps). Due to the introduction of 
the charges by means of the Poisson summation for- 
mula there are two imaginary-valued contributions to 
the Hamiltonian: the coupling of thermal charges n to 
quenched charges q and the coupling of vortices N to 
charges n. However, the partition sum is real- valued since 
the imaginary parts cancel in pairs exp[— Tl({N}, {n})} + 
exp[-H({N},{-n})}. 

Both gases are dilute for K 3> 1 and H -C 1, where 
the core energies 



E v 



1 
2H 



P 2 1 
AttK 



(14) 



are large (7 ~ 1.6 is a constant of order unityEl emerging 
from the lattice Green function) . Although these core en- 
ergies are determined by K and H we will consider them 
formally as independent parameters in the subsequent 
treatment. 

The partition function of this effective model (|l^) reads 
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Z = Z(2nK/p,E v ,E Cl {Q},{q}) = ^ e~ n . (15) 

{W},{n} 

The primed summation is constrained by neutrality con- 
ditions 



(16) 



R 



since the real part of the energy of nonneutral configura- 
tions diverges logarithmically with the system size. 



D. Replication 

The direct analysis of the partition sum (0) in the 
spirit of the KT renormalization group (RG) is prevented 
by the broken translation invariance in the presence of 
the quenched vortices and charges. Since we arc aiming 
at the calculation of disorder-averaged quantities, such 
as the correlation function (^), we may perform the dis- 
order average on the replicated partition sum already at 
this stage.LJ 

For an integer number of replicas n the Hamiltonian 
of the replicated system is most conveniently expressed 
in a vector notation, where n r = {n^} for a = l,...,n 
is a replica charge at site r and Nr = {-/V R } is a replica 
vortex at site R: 



Ti = Ti. v + H. c + Tt v 



T~tvc = 



l - NR-K v -N R ^ln|R-R'| 



R^R' 

]Tn r -e v -Nr 



R 

-y, 

2 V 



n r • K r • rv 2ir In I r — r' 



^n r • E c • n r , 



ip^w(R - r)N R • n r . 

R,r 



(17a) 
(17b) 

(17c) 
(17d) 



All parameters are incorporated in the replica-symmetric 
coupling matrices 



K, 



a b 



KS 

P 2 



ab 



a K 



1 + na T K ' 

L 



(J 



E? 

E ab 



J_, 
2H C 



ab 



ab 



ab 



(18a) 

(18b) 

(18c) 
(18d) 



Hamiltonian ( J17| ) is a generalization of that given in 
Sec. II of CO. In the limit a L — > oo we retrieve their 



(2.15). Furthermore, in this limit the contribution of the 
charge core energy proportional to a L enforces an addi- 
tional neutrality 



= 



(19) 



across the replicas, i.e., the weight of nonneutral 
configurations vanishes in the partition sum. Then 
the vector n can be decomposed into pairs of type 
(0, . . . , 0, 1, 0, . . . , —1, 0, ... 0) at every site. Such pairs 
correspond to the vector charges of CO. However, for 
finite a h the summation over charges is not restricted 
by this additional neutrality condition. Note that Eq. 
(|17| ) contains an explicit coupling between vortices and 
charges. 



E. Duality 

The-,|-e|ffective model under consideration is self- 
dual,E2ft9 i.e., it is invariant under a mutual exchange 
of vortices and charges with a simultaneous substitution 
of parameters. For the unrcplicated model ( |l^ ) these 
substitutions are 

2itK p 

R <-> r, N <r~* n, Q <-> iq, <-> - — — , E v <-» E r 



2wK 



(20) 



This duality implies a relation 

Z(27rK/p,E v ,E c ,{Q},{q})e 
Z(p/2nK,E c ,E v ,{iq},{-iQ}) 



+nj27TK/p,{q}) = 
,+H :c (p/27rK,{-iQ}) 



(21) 



for the partition sum ,113 which simplifies our subsequent 
calculations because it enables us to relate renormaliza- 
tion effects arising from vortex fluctuations to renormal- 
ization effects arising from charge fluctuations or vice 
versa. In this way PK obtained recursion relations for 
the bond-disordered model with field from those of RSN 
without field. We are going to use the same strategy 
starting from the corrected flow equations.Ea 

We will proceed within the replica formulation (|lj^ 
where the self-duality of the model is reflected by the 
substitutions 



R <-> r. N <-> n, K v <-> K c , E v <-> E c 



(22) 



model, i.e., Eq. (17c) becomes equivalent to their Eq. 



The excess contribution (12c) to the Hamiltonian of the 
unreplicated system effectively modifies the action of the 
longitudinal disorder distribution. The duality K v <-> K c 
implies a L ^ -(4n 2 K 2 /p 2 )a T /(l + no T K). The latter 
relation is consistent with Eq. (20) if one keeps in mind 
the presence of the excess term in Eq. (^l|) . This term 
gives rise to the denominator 1 + ncr T K that becomes 
unity in the limit n — *• that eventually will be taken. 
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III. RENORMALIZATION 

After the Coulomb gas representation of our model has 
been established in the previous Section we now turn to a 
scaling analysis for the relevance of vortices rand charges. 
We basically follow Kosterlitz' RG approachtil to the XY 
model and its recent generalization!!] to the XY model 
with random phase shifts. To account for magnetic fields 
we have to include charges and their coupling to the vor- 
tices. 

Following Ref. |2^ we introduce the notion of "types" 
of replica vortices and charges. We consider charges at 
site r and site r' as being of the same or of different type 
if n r = rv or n r ^ n r < in analogy to the vortex types. 
From the core energy contribution to the Hamiltonian 
( |l7| ) we identify fugacities 



y v (N) = e - N - E - N 

l, c (n) - e- n E * n 



(23a) 
(23b) 



for each type of particles. Because of the replica sym- 
metry of the matrices E v and E c different particle types 
have the same fugacity if they belong to the same "class," 
that is, if they are of the same type after a permutation 
of replicas and/or an inversion of the replica vector (such 
as n r — * — n r ). 



A. Scaling analysis 

In the limit of vanishing field all charge fugacities 
go to zero. The relevance of a weak field can be exam- 
ined by looking at the scaling behavior of the charge fu- 
gacities. As soon as the fugacity of some charge type 
increases under rescaling, the field constitutes a relevant 
perturbation. Analogously, the relevance of vortices sig- 
nals the importance of topologically nontrivial spin con- 
figurations. However, while a relevance of vortices indi- 
cates a tendency to reduced spin correlations on large 
scales, a relevance of charges expresses a tendency to en- 
hanced spin correlations. 

A rescaling of lengths R — > e~ rf/ R and r — > e~ dl r with 
an infinitesimal increase on a logarithmic length scale I 
has the following effect on Hamiltonian (0): the angles 
to are invariant but the logarithms generate additive con- 
stants, I n 1 . . . | — ► In | . . . | — dl. Using the neutrality con- 
ditions (|l6|) the twofold spatial summations over these 
constants can be reduced to single summations. Then 
the rescaling effectively amounts to an additive flow of 
the core matrices E, or equivalently to a multiplicative 
flow of the fugacities 



d 

li 

d_ 

Jl 



2/v(N) = (2 — 7rN • K v • N)j/ V (N), (24a) 



seal 



For a finite number n of replicas we can now determine 
the stability of QLRO to vortices and charges at given 
temperature K^ 1 = T/J and disorder strengths cr n by 
looking for the most relevant class of particles. We fo- 
cus here on the charge sector, since the vortex sector has 
been covered in Ref. ^2|. The most relevant charge class 
is the one that minimizes 



n K r 



P 

Air 2 



a \ a 



(25) 



The minimum is provided by charges of class j for 
<j L < K~ x and by charges of class fj, for a L > K~ x . 
The class "T" is composed of all replica charges n that 
contain an elementary charge in a single replica [e.g., 
n = (0, . . . , 0, ±1, 0, . . . , 0)] and the class "||" is com- 
posed of all replica charges n that contain an elementary 
charge in one replica and an opposite elementary charge 
in a second replica [i.e., the CO vector charges, such as 
n= (0,... ,0,4-1,0,. ..,0,-1,0,. ..,0)]. 

Let us consider for illustration vector charges that are 
composed only of s elementary charges, thereof Sf are 
positive and Sj are negative. The fugacity of such charges 
scales like 



y c (n) ~exp<l 



[(« T + s^K- 1 - (s T - Sl ) 2 a L ' 



2 1 



ill 



(26) 



y c (n) = (2 - nn ■ K c ■ n)y c (n). 



(24b) 



seal 



Figure H displays the stability boundary in the (K~ 1 ,cr L ) 
plane for various classes of charges, which are irrele- 
vant on the high-temperature side of the corresponding 
(full) lines. Charge class f is most relevant on the high- 
temperature side (<t l < K^ 1 ) and charge class fj. is most 
relevant on the low temperature side (<j l > K^ 1 ) of the 
dashed line. 

There is an important difference in the scaling analysis 
for charges and for vortices. The region of the phase di- 
agram, where all n-component charges are irrelevant, is 
independent of the number n of replicas. In addition, the 
boundary of this region is invariantly given for all integer 
n by the charge classes f or | J, (strictly speaking, the class 
t| exists only for n > 2). This independence on n indi- 
cates that the replica limit n — > can be taken in a very 
simple way. In contrast, the region where all vortices are 
irrelevant depends explicitly on n and so does the class of 
the most relevant vortices at lower temperatures (see Fig. 
2 of Ref. ^2|) . Therefore the replica limit will be easier to 
take for charges than for vortices. The save proceeding 
for both cases is to consider the collective renormaliza- 
tion of the physical parameters by all classes of vortices 
and charges, which will be carried out subsequently. 

We anticipate here that through this proceeding we 
obtain first closed expressions for the collective effects of 
charges, for which the perturbative expansion in small 
charge fugacities can be justified a posteriori. The vari- 
ous terms in this expansion can again be identified with 
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vector charge classes. The leading terms correspond to 
the most relevant classes identified above and the stabil- 
ity region (shaded area) of Fig. remains unchanged. 



B. Screening 

The renormalization of parameters on large scales due 
to fluctuations on small scales will now be calculated by 
integrating out pairs of particles fronp_the partition sum 
following the approach of Kosterlitz.El Since our model 
possesses the vortex-charge duality, it is sufficient to per- 
form the actual calculations only in one particle sector. 
The renormalization effects of the other sector are then 
obtained by a duality transformation. We adopt the re- 
sults for the vortex sector directly from Ref. How- 
ever, we extend these previous calculations at two points. 
First, we calculate the flow equations for the fugacities 
to higher order in the fugacities. This extension is neces- 
sary to consider the parameter renormalizations not only 
within the phase with QLRO, but also in the phase where 
magnetic fields are relevant. This requires flow equa- 
tions of higher order in the charge fugacities, which we 
can obtain through the duality mapping only from the 
terms of higher order in the vortex fugacities. Second, 
when the vortex dipoles are integrated out in the pres- 
ence of charges, they renormalize the interactions among 
the charges (and vice versa). This feature was naturally 
absent in Ref. ^2|, where fields were not included. 

For technical convenience we allow replica vortices to 
take positions in the continuous space and not only on 
the lattice. The lattice provides only a cutoff (taken 
to be unity) for the minimal distance between vortices. 
Our present goal is to integrate out all configurations 
in the partition sum, where the smallest distance be- 
tween two replica vortices (of any type) lies in the range 
1 < |R-i — R-2I < 1 + dL~-.We follow Young's treatment 
of a vector Coulomb gas,E2l which can be applied directly 
to our vectors N. During this procedure two cases have 
to be distinguished: (i) the two close vortices are of op- 
posite type, i.e., Ni + N2 = (in the continuous space 
subscripts to N denote the particle label instead of the 
position) and (ii) these two vortices are not of opposite 
type, i.e., N a + N 2 ^ 0. 

In case (i), where the pair of close vortices is composed 
of a vortex and its antivortex, they essentially annihi- 
late. The leading interaction with the other particles, 
which are approximately considered to be far away, is 
through the dipole moment of the pair. Because of the 
polarizability of this dipole moment the-»air screens the 
interaction between the other particles The screen- 
ing of the interaction among vortices can be captured by 
a flow of the coupling 



ab 



d 
Jl 



K v = 2tt 3 K v • C v • K v , 



(27) 



£ N a N b y 2 v (N) 



(28) 



is the replica vortex density correlation at the cutoff dis- 
tance [see Eq. (13) of Ref. |22j. Due to the angular cou- 
pling between vortices and charges there is an additional 
screening of the interaction among charges (see Appendix 



d_ 

dl 



(29) 



Besides the screening of these coupling matrices, the an- 
nihilation of such dipoles changes the free energy per unit 
area and per replica by 



d_ 

Jl 



7? ^ — ' 

scr N _^ 



2/v(N). 



(30) 



In case (ii) the two close vortices appear (in the eyes 
of the other particles, which are assumed to be far away) 
essentially as a single vortex of type N 3 = Ni + N 2 , i.e. 
the dominant interaction with the other particles is not 
given by the dipole moment of the pair but by the total 
vorticity. Therefore we can "replace" the two vortices 
by the single one, i.e., we transfer the statistical weight 
of the entire particle configuration to a different particle 
configuration, where the two close vortices 1 and 2 are 
replaced by a vortex of type 3 and all other particles are 
unchanged. By analogy with Young's analysisEII [leading 
to his Eq. (75)] we obtain 



d_ 

Jl 



2/v(N) 



Ni,N 2 #0 



^N,N!+N 2 yv(Ni)2/ v (N2). (31) 



(f) are preserved under such 



The neutrality conditions 
types of recombination. 



C. Flow equations 

In order to establish the RG flow equations, we com- 
bine the contributions from scaling and screening. Us- 
ing the duality relations ( p2| ) we infer the contributions 
by charge fluctuations immediately from those by vortex 
fluctuations that have been given above. Thus, after col- 
lecting all terms, we find the complete flow equations for 
the replicated system 



d 



(2-^N.K v .N)y v (N) 

+7T ^2 < ^N,N 1 +N 2 Z/v(Ni)y v (N 2 ), (32a) 

Ni,N 2 #0 



where 



(2 - 7rn • K c • n) y c (n) 

5n. ni +n 2 yc(ni)2/ c (n2), (32b) 



• E 

ni,n 2 ^0 



G 



and 



^K v = 2tt 3 K v C v K v - |p 2 C c , 
j t K c = 2tt 3 K c C c K c - |p 2 C v , 



(33a) 
(33b) 
(33c) 



n^O 



Since we establish the flow equations for an infinitesimal 
change of the cutoff, we assume that we can simply add 
vortex and charge contributions, i.e., that there are no 
terms mixing vortex and charge fugacities. Such a mix- 
ing could arise in principle if one implies a minimal cut- 
off distance not only among vortices and among charges, 
but also between vortices and charges. We are inclined to 
believe that a hard-core interaction between vortices and 
charges is not necessary, because we do not expect diver- 
gences related to short distances of that type: due to the 
imaginary- valued nature of angular coupling (f2d) in the 
Hamiltonian it seems plausible that configurations with 
close vortex-charge pairs give no essential contribution to 
the partition sum. 

Since not only the original couplings K but also the 
particle density correlations C are replica-symmetric, 
this symmetry holds also for the screened couplings. 
Even more, the renormalized matrices K v and K c can 
be reexpressed in terms of renormalized parameters K 
and cr 11 preserving relations (fl8^b,b): Equations (|33|a,b) 
are equivalent to 



dl' 



: K 1 = -2^ 3 C v , con + -p 2 K~ 2 C c , con , (34a) 
+a n K(2 + na n K)(C c , con + nC c , dis )]. (34b) 



Thereby the particle correlation matrices were decom- 
posed into the "connected" and "disconnected" contri- 
butions according to 

C ab = C con S ab + C dis (35) 

for both vortices and charges. 

D. Replica limit n — > 

So far, the derivation of the RG flow equations for 
the replicated system have been conceptually straight- 
forward. The less trivial part is to take the replica limit, 
i.e., to perform an analytic continuation for n — > 0. 

The first step in this direction was the reduction of the 
K matrix flow equation to the parameter flow equations 
(|34|). We focus on these equations, since the physical 
large-scale properties of our system show up in the flow 
of these parameters. In order to send n — > in these 
flow equations, we have to overcome only one hurdle: to 



make n an explicit parameter in the connected and dis- 
connected correlations. Here again we follow the route 
of Ref. |22|, where the problem was solved for the vortex 
correlations, and treat charges on the same footing. 

One easily realizes from the flow equations for the 
fugacities that large vortices (charges) with \N a \ > 1 
(\n a \ > 1) are less relevant than small vortices (charges) 
with |iV a | < 1 (|n a | < 1). Therefore it is a good approx- 
imation to restrict the further calculation to the small 
particles (this approximation is made for simplicity, the 
inclusion of large particles is possible). The correla- 
tions can then explicitly be evaluated after introducing 
a auxiliary random variable A. Its distribution is de- 
fined to be Gaussian with average [^4]_4 = and variance 
[A 2 ]a = 1/2 (averages over A are denoted by [. . We 
decompose the core energy matrices (|l8|c,d) into 



E 



ab 



ES 



ab 



E 



(36) 



for vortices and charges (one specifies by adding the cor- 
responding subscripts "v" or "c" ) . The initial values of 
E given in Eq. (|l4|) and 



EL = 



ir^a K 2 
1 + n<r T K ' 



E r = - 



7P 2 cr L 

47T 



(37) 



Next, we introduce weights (proportional to the fugaci- 
ties e~ 2E of the elementary particles) 



z± 



-2(E±AV2E) 



1 + Z + + Z- , 



(38) 



and, combining Eqs. (|23|) with (|28|) and its dual coun- 
terpart, we rewrite the correlations as 



C ab 



[(z+ + z-)z 

[[(* 



z_) 2 z"- 2 l 



(39a) 
(39b) 



where we take a ^ b. In this representation the number 
of replicas became an explicit variable that can be sent 
to zero, resulting in 



Cdis — — 
Ccon — 



24 



1 + Z+ + Z 

z + + Z- + 4z + z_ 



J .4 



(l + z + + zSf 



(40a) 
(40b) 



Then the replica limit of the flow equations (|34|) for K 
and <r n can be taken right away, since now all depen- 
dences on n are explicit. 

The contributiong-.to the flow of the free energy (|33|c) 
can be rewritten asE3 



-Y / yl(^) = -K~i]A^[inz v } A , 



(41) 



N#0 



where in the last expression the replica limit has already 
been performed. The analogous charge contribution is 
obtained by duality substitutions. 
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E. Fugacity expansion 

Combining expressions (§|c), @, ©, and (||) the 
flow of the physical quantities K, a n , and J 7 can be ex- 
pressed through the fugacities z± for both types of par- 
ticles. According to Eq. ( [Til) these fugacities are small 
for large K and small H and one might be tempted to 
simplify the expressions (40) and ( f4l|) by a truncated ex- 
pansion in these fugacities. _____ 

The early studies of the disordered XY modcMI___d 
made use of this fugacity expansion. The invalidity of 
this expansion was recognized later on and was consid- 
ered as a possible reason for the destj__|Clion of QLRO 
even by infinitesimaLta-. W*ak disorder However, it 
was shown recently-*-.—! for the model without fields 
that the mere breakdown of the fugacity expansion does 
not imply the relevance of vortices and the destruction 
of QLRO. In other terms: in the vortex sector the closed 
expressions (^0|) and ( ^l| ) are finite and they are to be 
evaluated nonperturbatively. Therefore it is a central is- 
sue of the present work to reexamine the validity of the 
fugacity expansion for charges that underlies Refs. [l6 18 

To be explicit let us examine the "smallness" of the 
charge fugacities z c ± defined in Eq. (|38|). In contrast to 
the real- valued vortex fugacities they are complex- valued 
since \f~E~ c is purely imaginary. Their absolute value 



z c ± 



-2E C 



-1/H 



Hi 



(42) 



is small for weak fields independently of the disorder rep- 
resented by the variable A. Therefore it is safe to ex- 
pand the expressions correlations ( |40| ) for charges before 
performing the average over A. Again, this is in sharp 
contrast to the case of vortices, where the fugacities be- 
come arbitrarily large in the tails of the distribution of 
A. Thus, the mechanism that led to the failure of the fu- 
gacity expansion for vortices is not effective for charges. 

Thus, we are entitled to use the fugacity expansion for 
charges to study the stability QLRO to weak fields. We 
expand the charge correlations ( |io| ) [and analogously the 
source of the flow of the free energy, the dual counter- 
part of Eq. ((4l|)] to second order in the fugacities z c ± 
and perform the average over A: 



C c ,dii 
Cc_c.cn 



= -2y, 



4j/cTT 



[lnz c ] A = 2y 2 c] -y 2 cU - y 2 ^ 



(43a) 
(43b) 
(43c) 



We have replaced the terms in the series by vector charge 
fugacities using the identification ( p3| ) and substituted 
the symbols y c (n) by the more illustrative arrow sym- 
bols. The initial values of these fugacities are given by 



All bare charge fugacities are determined by only three 
physical parameters K , a L , and H . This relation is aban- 
doned under renormalization, where the flow of K and 
a L follows Eq. (|34| ) and where we have a whole set of 
flow equations fl32|) for the charge fugacities instead of a 
single flow equation for H . It is interesting to note that 
a fugacity with s arrows is of order Hq in the original 
field [remember relation (pj])]. This order s = + Sj 
is just the total number of positive /negative elementary 
charges introduced in Eq. (Eq) above. 

Since the fugacity expansion is valid for the charges, 
it is legitimate to analyze the stability of QLRO by de- 
manding that all classes of vector charges have to be ir- 



[II A 



relevant. This analysis has been performed in Sec 
leading to the shaded stability area in Fig. |. As long 
as all charges are irrelevant, the second order terms in 
the flow equation (|3-jb) play actually no role. For a qual- 
itative study of the flow of charge fugacities slightly in- 
side the region where charges are relevant, we keep only 
the two charge classes y c | and y c u- This should be a 
reasonable approximation since it is always one of these 
two classes that brings about the instability of QLRO at 
the low-temperature border of the stability region and 
that represent the most relevant charge class among all 
classes and for all parameters, as discussed after Eq. (|26|). 
Therefore we truncate the charge fugacity flow equations 
_ 3) to 



y cT = [2 - (p 2 /4^)(K- 1 + a L )]y c] - 2^ cT2Mi , (45a) 



dl 



dl 



y cll = [2 - (p 2 /27T)K- 1 ]y cU + 2n[y 2 ^ - 2y 2 ctl }. (45b) 



However, we stress that these truncated flow equations 
are incomplete, i.e., these two selected fugacities are cou- 
pled to all other fugacities in the complete flow equations. 
We have dropped also charge class |f , which also brings 
about only quantitative modifications. Only for a L = 
they are as relevant as charges fj an d should be kept, 
e.g., to obtain the vanishing of the disconnected correla- 
tion (g|a). 

Concerning the processing of the vortex sector, we ac- 
tually don't go beyond Ref. [22[ That is, we restrict our 
quantitative analysis to the region, where vortices are 
irrelevant and the second-order terms to the vortex fu- 
gacity flow equations (|3^a) can be discarded. The linear 
flow equations for all y*.(N) can be reduced to flow equa- 
tions for E v and E V E3 Avoiding a fugacity expansion, 
the vortex correlations C v can be evaluated asymptoti- 
cally for large length scales and can be expressed in terms 
of an effective vortex fugacity y v , 



2/cT 
2/cU 



-[l/2H+(-yp 2 /47T){K- 1 +cr L )] 
-{1/H+{~ 1P 2 /2tt)K- 1 } 



„ - p -[l/H+(-fP 2 /^)(K- 1 +2a L )] 

i> c TT — e 



(44a) 
(44b) 
(44c) 



C — — 9t*ii 2 

^v.con ^ > y v i 



C v 1 



-2(1 



[In; 



2 

7* 



(46a) 
(46b) 

(46c) 
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with a fraction of "polarizable dipoles" 

t* =mm(l,l/2a T K). (47) 

This fraction is unity for temperatures K^ 1 > 2a T and 
goes to zero as the temperature drops below K^ 1 < 2er T , 
where vortices start to freeze. The initial value and t 
renormalization flow of the effective fugacity is given b; 



IV. RESULTS 



1 



z + + Z- 



Z + + Z- 



dl 



y w = [2--KT*K{l-a i r*K)]y v . 



(48a) 



(48b) 



To summarize the above derivations we collect the full 
set of flow equations that constitute one of our main find- 
ings: 

^ If -1 A 3 * 2 2 j^-2 2 

—K = 4tt r y v - irp K y cV 

-K- 2 y 2 cU ], 



-y v = [2 - t:t*K(1 - a T T*K)]y v , 



(49a) 

(49b) 
(49c) 



-y c] = [2 - (p 2 /4 7 r)(^- 1 + a L )}i M - 27ry cT2/cU , (49d) 



-y cU = [2 - {p 2 /2n)K- 1 ]y cU + 2^ - 2y c 2 u ], (49e) 

(49f) 



^ = --^i/v-n 2 »cT-Vcu) 



This set of flow equations has no longer an explicit du- 
ality between vortices and charges (even if we ignore the 
second-order terms in the flow equations of the charge 
fugacities), because the approximations used to treat the 
coupling of vortices and charges to disorder have to be 
taken in a fundamentally different way, mainly because 
of the validity (invalidity) of the fugacity expansion for 
charges (vortices). In particular it was possible to retain 
only two charge classes, whereas the collective effect of 
all vortex classes is represented by the effective vortex 
fugacity and the appearance of the parameter r* (that 
actually parametrizes the most relevant vortex class). 

However, if we specialize flow equations ( ff9| ) to the 
case without disorder, the duality reappears and we re- 
trieve the flow equations of Ref. p"3p3 If we switch off 
the field, i.e., if we drop the charge sector, the flow equa- 
tions (^9|) are re duce d per construction exactly to those 
derived recently.! 2 tI 23 ! In the presence of random fields 
[where a L = oo and y c j = 0, see Eq. (fti})], the flow, 
equations reduce in the charge sector to those of CO.t3 
However, the present flow equations allow also the con- 
sideration of finite a L . The charge sector for uniform 
fields has structural similarities with the flow equations 
that havejpeen given (but that have not been evaluated) 
in Ref. ||.Ej The flow equations in Ref. [l8] miss the contri- 
butions from fugacities y c \\ and the freezing of vortices 
for t* < I. 



We now determine the large-scale properties of the 
model by integrating the flow equations (|49|). Thereby 
the initial parameters K and cr 11 are mapped for I = oo 
onto Koo and c^. The fugacities, which are initially 



given by Eqs. (44) and (|48ja) , flow independently. Note 
that t* is also a flowing parameter which is determined 
by the flowing K and a T through relation (0). 

From the flow equations one easily reads off how vortex 
and charge fluctuations renormalize the physical param- 
eters qualitatively. Equation (^a) means that the effec- 
tive temperature is increased by vortices but reduced by 
charges, corresponding to disordering tendency of vor- 
tices and the ordering tendency of charges. Vortices al- 
ways increase the strength of disorder but the effect of 
charges is ambivalent in Eq. (|49"|b) . The competition be- 
tween the various terms leads to a rich structure of the 
phase diagram. For the clarity of the analysis we first 
break the problem into parts that we eventually assem- 
ble to the complete picture. 



A. Vortex matter 



For completeness and further reference we recall 
the main . .properties of the system in the absence of 
fields :E3c3'E^I At temperatures below the KT transition 
temperature of the pure system, that is, K^ 1 < tt/2, 
vortices are irrelevant even in the presence of random 
phase shifts with a strength below the critical value (see 
Fig. § 



i 



^ for | < K-J < |, 
for K^<\. 



(50) 



The increase of K^ 1 and a T under renormalization im- 
plies that the phase with QLRO is less extended in the 
(K~ 1 ,a T ) plane of the unrenormalized parameters as 
compared to the (K^ 1 ,cr^,) plane of the renormalized 
parameters. Nevertheless, this shrinking is only a quan- 
titative modification and QLRO is stable for low temper- 
atures and weak disorder (see Fig. 3 of Ref. p2| ). 

It is evident from the dependence of the flow equations 
on the parameter r* that the physics within the phase 
with QLRO is different for high temperatures and for low 
temperatures. Vortices start to freeze and the breakdown 
of the fugacity expansion sets in at temperatures below 
the line K^. 1 = 2aJ a separating these two regimes. 

The spin correlation function T p decays algebraically 
as in the spin wave approximation (|^) but with a renor- 
malized exponent 



p-_ 

2tt 



K 



(51) 







Through the renormalization of K and a L the value of 
r\oo depends implicitly also on the transverse disorder. 
Its value is r\oo = p 2 /4 at the high-temperature end 
of the transition line and drops continuously down to 
Voo = p 2 /16 at the zero-temperature end of the transi- 
tion line. The exponent 



Vc 



1 



(52) 



of the Edwards-Anderson spin correlation (j^) assumes 
values in the range < rj < 2 and also varies nonuniver- 
sally along the QLRO/SRO transition. When the tran- 
sition line is approached from the high-temperature side, 
the correlation length diverges exponentially (with excep- 
tion of a special point at K^ 1 = 2cr^E3) with decreas- 
ing distancp-ifcom the transition line as in the absence 
of disorder .E3I23 Within the QLRO phase glassy features 
are found in *»rtex density correlations for temperatures 



one might question again whether charge classes different 
from the two considered ones could change the value of 
x. This is not the case and can be understood from an 
inspection of the x values for the other classes. These val- 
ues are obtained by dividing the exponent in Eq. ( p6| ) by 
the power s — s-\ + Sj of the scaling relation between the 
fugacity and the field. The different order of the charge 
fugacities in the field is also the reason that the separa- 
tion line K^ 1 = a L , where dominance switches between 
charges of class j or tj., does not coincide with the sep- 
aration line a L = Air/p 2 where the exponent x has its 
nonanalytic dependence on K~ x and <r L . 

From the assumption that the magnetization yields 
an extensive contribution to the free energy, AF ~ 
L 2 ~ Hq^ x and the definition of the magnetization 
M ~ OAF/ 8Hq we identify the exponent 



s- -1 



for at<^, 



for a^> 



in 



(56) 



B. Charge matter 

Let us now consider the complementary case where 
we evaluate the effect of fields in the absence of vor- 
tices. Taking into account renormalization effects, the 
linearized charge-fugacity flow equations indicate the ir- 
relevance of fields at high enough temperature even for 
arbitrarily strong disorder: 



K~ 



> 



^-cr^ for cr£ < %- 



1 7; 
P 2 



for 



> 



4tt 



(53) 



In the low-temperature region of the phase diagram, 
where the magnetic field is relevant, we can determine 
the critical exponent S that defines the scaling relation 
M ~ H^ S between the magnetization M and the field 
Hq. From the linearized flow equations of the charge 
fugacities we determine the exponent x in the scaling 
relation Ho ~ L x between the field amplitude and the 
system size L = e l . From the two charge classes under 
consideration wc find 



j2-{p 1 l^)(K- 1 +a L ) 



Hn 



Hi - j L 2 -(P 2 /27r)i<"-\ 



(54a) 
(54b) 



From these two scaling relations one obtains at first sight 
two different values for x. The actual value of x, which 
describes the relevance of Hq, is given by the larger of 
the two values. Therefore we find 



2 -!?(*- 



1 



l-K- 



<J L ) for <J L <%, 
for a L > p, 



(55) 



as special cases of Eq. (26). To determine this expo- 
nent it is important to keep in mind that different vector 
charges are of different order in the field. At this point 



This expression reproduces the universal value 5 = 15 
for fields, with p = 1 at the KT transition of the pure 
system.U 

At low temperatures, where the field is relevant, we 
find two different phases. For weak disorder a L < 4ir/p 2 
the ordering tendency of the fields dominates over the 
disorder and LRO is established. In the opposite case 
<j L > 4ir/p 2 the disorder wins over the field and results 
in quasi-short-range order (QSRO) characterized below. 
The phase diagram, obtained from a numerical integra- 
tion of the flow equations in the absence of vortices, is 
depicted in Fig. |]. 

LRO is found in a region a < An/p 2 and Kq 1 < 
8tt/p 2 — <j l of the unrenormalized parameters. Under 
renormalization both charge fugacities saturate at a fi- 
nite value and drive these parameters to zero, 



2/cT 



1 

7T 
1 

! 

7T 

n 



K « —v 
p z l 

~pH- 



(57a) 
(57b) 
(57c) 
(57d) 



On large scales K~ Y and a L decay only inversely propor- 
tional to the logarithm of the length scale, and so does 
the effective exponent r\ « l/l. Therefore this fixed point 
corresponds indeed to a phase with LRO, where 



r P (r) 



const. > 



(58) 



for r — > oo. The fact that the fugacities saturate at a fi- 
nite value 1 j~K is consistent with a very (maybe infinitely) 
strong field Hq at the fixed point, since even an infinitely 
strong initial field would lead to finite fugacities, see Eq. 
(E3). The analysis of the LRO phase cannot be taken 
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quantitatively, since the flow of K^ 1 and <r L entails the 
relevance of other charge classes that have been ignored 
in this analysis. However, we expect modifications to be 
only quantitative since the omitted charges are less rel- 
evant than the retained charges [see the discussion near 
Eq. ©]. 

QSRO is found in a region a > An/p 2 and Kq 1 < 
Air/p 2 of the unrenormalized parameters and identified 
with the following asymptotic flow: 



?M -> 0, 
2/cU 



P 



2ir \ AmK n 



K' 1 -> K~} < 



P 



AixK 2 



47T 



1 



4nK n 



I. 



(59a) 
(59b) 

(59c) 
(59d) 



Since both <j n increase without bounds, there is no 
fixed point i n a s trict sense. According to our anal- 
ysis in Sec. 



Ill A 



all charge classes different from |J, 
become actually irrelevant at temperature not too far 
below K^ 1 — An/p 2 . However, at significantly lower 
temperatures other classes will be relevant, e.g., ttll at 
K^. 1 — 2tt/p 2 . The divergence of a L means that the sys- 
tem behaves asymptotically like in the presence of ran- 
dom fields. Since both a 11 don't couple back into the flow 
equations, their effect on the angle difference correlation 
function can be evaluated in the spin-wave approxima- 
tion accounting for the scale dependence of a L , 



{(Or 



1 



2 V 2ttK 2 



P 



1 



P 



AnKl 



ln 2 r, (60) 



as found in Refs. ^6|3^j36j. In the context of crystals 
grown on disordered substrates this behavior of the cor- 
relation function is called "super-rough" .□ To the ex- 
tent that only the class || contributes to the renormal- 
ization, K is not renormalized. For the random field 
model this non-renormalization can be shown, in gen- 
eral on the basis of a statistical symmetryEj'LH Then 
the amplitude of the difference correlation in front of 
In 2 r is universal in terms of the transition tempera- 
ture Tjj, which is determined by K^ 1 — Air/p 2 , namely 



,) 2 ) » (2/p 2 )(l - T/T U ) 2 In 2 r for T < T n (Ox 



quantitative analysis has to be restricted to parameters, 
where the fixed-point value of y^ is small. Otherwise 
higher orders in this fugacity can no longer be neglected). 
Our calculation agrees in the pxefactor of In 2 r with that 
of Carpentier and Le Doussal.Lj The increase of the dif- 
ference correlation faster than logarithmically but slower 
than linearly in the length implies a faster than algebraic 
but slower than exponential decay of the spin correlation, 



r P (r) 



,-(l-T/T c ) 2 lnr 



(61) 



i.e. the correlation length is still infinitely large. To em- 
phasize this we call the order quasi-short-ranged. 



In the QSRO phase the competition between the spin 
coupling energy and the magnetization energy leads to a 
state, which has less order than each of the contributions 
alone would have in the absence of the other contribu- 
tion (QLRO or LRO, respectively). This competition 
results in a frustration typical for glassy disordered sys- 
tems. Since the symmetry-breaking field does not induce 
a finite magnetization, the scaling relation and the expo- 
nent 5 are actually meaningless in this phase. 

The QSRO phase can be considered as "glass" since 
the disorder (a L ) dominates the physics on large scales. 
In addition, the dominant fluctuations are vector charges 
of type tl, i-e. fluctuations that are correlated in differ- 
ent replicas and usually would imply a finite Edwards- 
Anderson order parameter. However, we are not aware of 
such an order parameter in the strict sense for this phase. 
But we observe that the Edwards-Anderson correlation 
(||) decays only algebraically with finite rj^ , whereas the 
spin correlation T p decays faster than algebraically. This 
indicates that the conformations of 9 are disorder domi- 
nated and that thermal fluctuations are not stronger than 
in the absence of disorder. Additional evidence for the 
glassy nature of this phase can be derived from dynamical 
properties£j which are beyond the scope of the present 
work. 

The transitions between the three phases LRO, 
QLRO and QSRO are of different nature. When the 
LRO/QLRO transition is approached from above, the 
spin corcelation exponent reaches the universal value 
r^oo = 4,L3 and <5 = when it is reached from below 
since the field becomes irrelevant. These values are thus 
independent of weak disorder with < Air/p 2 . Along 
the QSRO/QLRO transition the exponent 7700 > 4 has 
a nonuniversal dependence on the disorder strength, and 
5 is not a meaningful quantity However, rj^ = A/p 2 is 
universal along this line £3 

The phase transition LRO/QSRO is located near er^ = 
Air/p 2 , see Fig ||. The transition line found from the nu- 
merical integration of the flow equations is bent towards 
smaller a L and the bending decreases with decreasing 
field amplitude. In Eq. (|5^) we have found the scaling ex- 
ponent of the field to change along the line — Air/p 2 , 
which indicated that the instability of QLRO against the 
field differs above and below this line. However, a change 
of the exponent for initial scaling of the field does not 
necessarily imply that the renormalization flow drifts to 
different fixed points. We can provide an additional ar- 
gument in support of this location of the transition line. 
Let us identify the difference between the phases with 
LRO or QSRO by the presence or absence of y c | on large 
scales. We further consider very weak fields such that K 
and a L are very weakly renormalized (on intermediate 
scales). Then we may use an adiabatic approximation 
for the charge fugacities, i.e., we suppose the values of 
the fugacities to be determined through dy c ^/dl = and 
dy c fi/dl — as a function of slowly varying K^ 1 and 
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cr L . From the flow equations 
ent solutions 



one obtains two differ- 



LRO 



QSRO 



2/cU ~ 2¥ 



y c 2 T 



1 

yen 



p 

4tt 



p a 



Veil 
Vc\ = 



J_ 

2tt 



47T 



JttA 



0. 



(62a) 
(62b) 



which are precursors of the LRO and QSRO fixed points. 
Eventually these precursors flow towards the true fixed 
points because of the scale dependence of K and a L . 
These two solutions merge again at the line a L — 4ir/p 2 . 
From a linear stability analysis of the fugacity flow equa- 
tions one finds that the LRO precursor is stable and the 
QSRO precursor is unstable for cr L < Air/p 2 and the 
QSRO precursor is stable for a L > 4n/p 2 . The LRO pre- 
cursor is unphysical for a L > Air/p 2 since there y% < 0. 
Although this analysis supports the horizontal location 
of the LRO/QSRO phase boundary in the limit of weak 
fields, we should like to remind again that the transi- 
tion can be described only qualitatively since this analy- 
sis includes only two charge types and the other neglected 
types are expected to give quantitative corrections even 
to the borders of the LRO phase, where the fixed point 
values of the charge fugacities remain finite. 

It is worthwhile to point out that the irrelevance of 
fields in the present context is equivalent to the irrele- 
vance of the periodic crystalline potential for the rough- 
ening transition of crystals with correlated substrate dis- 
order (the mapping be twee n these two problems is the 
one discussed in Sec. [IB, <f) corresponds to the sub- 



strate conformation). The latter problem has been ex- 
amined in Ref. ^ and the phase diagram obtained there 
(Fig. 2 therein) is the analog to the onset of relevance 
of fields in our Fig. |2[ In that work the longitudinal 
bond disorder originated from correlated <f> as in Eq. (|j]) . 
Therefore the LRO in the present context actually corre- 
sponds to a logarithmically rough surface profile that is 
locked into a configuration parallel to the disordered sub- 
strate. The present LRO/QSRO transition corresponds 
there to a transition from a locked-in and logarithmically 
rough surface to a super-rough surface. 

Recently Horovitz and Golub have examined the 
same problem, in view of its implications for Joseph- 
son junctions^ Using a Gaussian variational replica ap- 
proach they found a topologically similar phase dia- 
gram. But there are the following differences: The 
phase boundaries do not coincide exactly, since their ap- 
proach does not take into account parameter renormal- 
izations. In the phase where we find QLRO and LRO 
they also have QLRO and LRO. However, their phase 
with LRO is split into two sub-phases, their "Joseph- 
son" phase (where charges fj, are irrelevant) and "Coex- 
istence" phase (where charges fj. are relevant). We find 
charges fj. to be relevant in the whole LRO phase, i.e. the 
whole phase with LRO has the character of the "Coex- 
istence" phase. The boundary between their sub-phases 



coincides with the line K _1 = a L , where charges of type 
t| become relevant according to the linear scaling anal- 
ysis. However, both charges are relevant in the whole 
LRO phase, their appearance is triggered by the second- 
order contribution to the flow equation (f49|f) of y c ]i- 
The most important difference is the region where we 
find QSRO and they still find QLRO that qualitatively 
influences the dependence of the critical Josephson cur- 
rent on the junction area. Again, this difference is due 
to the presence of second-order terms in the flow equa- 
tions, now terms of order y^, to the flow equation ([49|f). 
The disagreements we find here are well-known from the 
random field XY model and have to be ascribed to the 
insufficiency of the Gaussian variational ansatz to cap- 
ture the physics represented by second-order terms to 
the flow equations. 



C. Complete picture 

We now discuss the phase diagram in the presence of 
vortices and charges. We distinguish the cases of uniform 
fields and random fields, since their effective longitudinal 
bond disorder is crucially different. 



1. Uniform fields 



T — 



For definiteness we consider the case a L = a 
of random phase shifts that are uncorrelated on different 
bonds. The general case with finite a L ^ a T deviates 
only quantitatively. As in the absence of disorderj!3 a 
phase with QLRO is stable only for p > 4. For weak 
disorder there is a LRO phase at the low-temperature 
side of the QLRO phase. For 4 < p < 4^/2 « 5.6 LRO 
is stable up to = tt/8, whereas for p > A\/2 LRO 
is stable only up to — 4ir/p 2 . Figure || depicts the 
phase diagram for p = 8, which is generic for all inte- 
ger p > 4\/2. As compared to the absence of vortices, 
one main effect of vortices is the restriction of the QLRO 
phase at higher temperatures, where SRO sets in. A 
second important effect is that the phase with QSRO is 
turned into a phase with SRO. Since a increases loga- 
rithmically with the length scale in the absence of vor- 
tices, it reaches the value a = tt/8, where vortices be- 
come relevant. On small scales vortices practically do 
not contribute to a renormalization of the parameters. 
If we start from bare parameter values at a temperature 
slightly below the QSRO/QLRO transition, only charges 
tl are initially relevant and increase a according to Eq. 
( p9| ) . Then vortices become relevant only on an exponen- 
tially large scaled 



L v w exp 



(1-T/T U ) 2 



(63) 



On this length scale SRO sets in, that is, L v repre- 
sents the correlation length of the system. At disorder 
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strengths a < 7r/8 the phase diagram, judged by its order 
at largest scales, appears to be reentrant, i.e., a phase 
sequence SRO/QLRO/SRO can be obtained by a pure 
temperature change. However, the difference between the 
low-temperature and the high-temperature part can be 
substantial in terms of the correlation length and prob- 
ably also in terms of dynamic properties. In a remote 
analogy the high-temperature part compared to the low- 
temperature part of the SRO phase should be as different 
as a gas compared to a very viscous liquid. 

Due to the restriction of our analysis to small fugacities 
of vortices and charges we can not rule out the possibil- 
ity of an actual transition between two different phases at 
a oo ~ for if" 1 < 4ir/p 2 . Since the boundary of the 
QLRO phase, which is well described in our approach, 



has a discontinuous slope at a c 



for K- 



4ir/p 2 



one might expect that this point is a tricritical point that 
should be merged by and a third transition line. 

For p < 4\/2 the phase boundary LRO/SRO is deter- 
mined by the competition of vortices and charges, which 
are both relevant already on smallest scales. The location 
of this transition can therefore not be determined. 

The phase diagram differs from the one found in 
Ref. [18] crucially at low temperatures, where the 
SRO and QLRO phases reached down to the point 
T = and a — and a reentrant phase sequence 
SRO/QLRO/LRO/QLRO/SRO was possible by chang- 
ing temperature. Although we still find a reentrancc 
SRO/QLRO/SRO the QLRO phase (as well as the LRO 
phase) has now a convex shape in the {K^, a) plane. The 
discrepancy can be traced back to the overestimation of 
vortex fluctuations in a vortex fugacity expansion, which 
is avoided in the present work. 

The phase diagram Fig. |^ has been obtained by a nu- 
merical integration of the flow equations (49) for two dif- 
ferent amplitudes Hq = 0.1 and 0.01 of the field and for a 
fugacity parameter j = 1.6 appropriate for the mapping 
of the original XY model onto the Coulomb gases. We 
find only a slight dependence of the phase boundaries on 
weak fields Ho <C 1. One could in principle evaluate the 
flow equations also for larger fields and would probably 
find a direct crossover from the glassy low-temperature 
SRO region to the nonglassy high-temperature region as 
suggested by CO. However, we do not pursue this issue 
further since for Hq > 1 it is certainly no longer legiti- 
mate to neglect the less relevant charge classes. 



2. Random fields 



The mapping discussed in Section II B allows us to 
study random fields with spatially uncorrelated <p by con- 
sidering a L = oo in the flow equations. In this limit the 
contributions proportional to y c -\ simply disappear from 
the flow equations, in which a L then no longer enters and 
the flow of a L can be identified with the flow of a L after 
subtracting the transformation term A/p 2 . 



Randomness in the field induces only a few changes in 
comparison to the case of uniform fields. The schematic 
phase diagram for the system with random fields and ran- 
dom phase shifts is depicted in Fig. |l]b, which is deformed 
only slightly by finite vortex fugacities and field ampli- 
tudes. In contrast to the case of uniform fields, the LRO 
phase disappears and the low-temperature QLRO/SRO 
boundary is located at K^ 1 = Air/p 2 down to a_= 0. 
Thus a phase with QLRO now exists for p > 2y/2H The 
spin correlation function has a finite exponent ( |5l|) (in- 
volving the original a L and not a L since we are interested 
in the correlation of 9 and not of 6). 

For temperatures T < Tj| the Cardy-Ostlund charges 
t| become relevant and still generate transverse bond dis- 
order according to Eq. ( |59"|) (which also describes the fur- 
ther decrease of longitudinal bond disorder). Therefore 
vortices become relevant again on scales (E33l) for T < Tj j_ . 



V. DISCUSSION 



To summarize, we have performed a renormalization 
group analysis of the model (|l|) in the Villain approxi- 
mation that allowed us to represent the model by two 
coupled Coulomb gases. The breakdown of the fugac- 
ity expansion for the vortex gasOnJ was treated nonper- 
turbatively (the flow equations contain the effective vor- 
tex fugacities y v that are nonperturbajliwe--fpictions of 
the original vortex fugacity ~ e -1 /^) «£3e3'c3 We have 
shown that an analogous breakdown does not occur for 
an expansion in the charge fugacities and have identified 
two most relevant classes of vector charges out of an in- 
finite set. The renormalization group flow equations for 
the interacting Coulomb gases have been established and 
evaluated for the limit where both gases are dilute, i.e., 
for low temperatures and weak fields. 

From these flow equations we have determined the 
structure of the phase diagrams, in particular the do- 
mains of stability of the phase with QLRO. Due to the 
breakdown of the fugacity expansion for the vortices this 
phase was found jtoJae more stable against disorder than 



found previously£3'E3 The critical exponents 7700 , rj^ , and 
5 have been given in Eqs. (|5|), (|5^), and (|5q) . 

In a strict sense our approach is valid only for the char- 
acterization of the QLRO phase, where vortices and fields 
are asymptotically irrelevant. Nevertheless we have at- 
tempted to complete the picture of the phase diagram by 
examining how vortices and fields become relevant, which 
was taken as indication for the LRO, QSRO, or SRO na- 
ture of the neighboring phases. Thereby we were also 
able to locate qualitatively the LRO/QLRO transition in 
the case of uniform fields. This transition was found to 
merge one corner of the phase with QLRO. Its location 
is in qualitative agreement with that found in a recent 
Gaussian variational approach^ where vortices had not 
been considered. On the high-temperature side of the 
QLRO phase vortices become relevant on relatively short 
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length scales and immediately lead to an increase of the 
effective temperature and a suppression of the charge fu- 
gacities. On the low-temperature side of the QLRO phase 
charges become relevant first. For weak longitudinal dis- 
order (<t l < An/p 2 ) the flow converges to a LRO fixed 
point, otherwise the flow tends on intermediate scales to- 
wards the QSRO fixed point before it eventually crosses 
over to the SRO fixed point. Since the crossover from 
QSRO to SRO occurs at finite values of the charge fu- 
gacities, this behavior might be modified by higher orders 
in the charge fugacities, which we are not able to incor- 
porate systematically at present. This difficulty excludes 
also the calculation of phase diagrams for uniform fields 
with p < 4 and random fields with p < 2y2 in the present 
framework. 

CO have suggested the possibility of a first-order 
transition between the glassy low-temperature SRO and 
the based on the nonglassy high-temperature SRO. To 
our understanding this scenario was based on low- 
temperature properties of the vortex fugacity flow equa- 
tion that are absent in the present analysis, where the 
vortex fugacity expansion has been avoided. 

We speculate that there could be a transition from a 
low-temperature SRO phase to a high-temperature SRO 
phase that merges the upper left corner of the QLRO 
phase. It could be possible that this transition line per- 
sists even at very strong disorder and ends up at a = oo 
and T = 0. This point corresponds to the so-called gauge 
glass model that represents a zero-temperature critical 
pointo However, for such a transition line we have no 
evidence other than the fact that corners usually occur 
only at a three-phase coexistence and that the gauge glass 
fixed point should have a continuation for weaker disor- 
der (finite a). It would be interesting to examine this 
scenario further, which probably can be achieved only by 
numerical methods. 

The studied model has physical applications in various 
fields. The most direct application is to XY magnets, 
where random phase disorder can be provided by non- 
magnetic impurities through the Dzyaloshinskii-Moriya 
interaction and-Jocal crystal fields break the spin rota- 
tion symmetry.EJ 

In the case when vortices are excluded, the model 
also describes crystal surfaces grown on a disordered 
substrate.QEa Such disorder can also be related rto dis- 
locations that leave the crystal at the surface.c2l The 
LRO/QSRO transition described above corresponds to a 
disorder driven transition from a rough crystal surface to 
a super-rough phase. In the context of planar Josephson 
junctionstl LRO corresponds to a finite critical current 
density even for arbitrarily large contacts and QSRO cor- 
responds to a critical current density that decays faster 
than algebraically with the contact area. 

The phase diagrams obtained here should also be qual- 
itatively similar to those_nf crystalline films on periodic 
substrates with disorderEru and of vortex latticesH The 
transitions occur between phases that differ in the degree 
of structural order. However, these systems are more 



complicate to analyze in a similar framework since they 
have a two-component displacement field instead of the 
one-component phase of model (|lj). In addition, the pe- 
riodicity of the lattice implies that both the particle in- 
teractions and the pinning energy have the same wave- 
length, i.e., p = 1. Uniform and disordered fields then 
correspond to pinning that is commensurate or incom- 
mensurate with the crystal. Since p = 1 it appears likely 
that there is no phase where the effects of pinning and 
the proliferation of dislocations are irrelevant on large 
scales and that the strongly relevant fields and disloca- 
tions cannot be well described by a direct generalization 
of the present framework. 
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APPENDIX A: THE COULOMB GAS MAPPING 

In this appendix some intermediate steps in the deriva- 
tion of the Coulomb gas model ( |l2"| ) aregiven. The start- 
ing point is the effective Hamiltonian pCf ) . The first step 
is to replace the scalar field m by the charge field n with 
the use of the Poisson summation formula: 

H = {y ( V(9r - A r + 27rM r ) 2 

r 

+ ^n 2 r +ipn r r y (Al) 

In the corresponding partition function all € [0, 2ir] are 
integrated over and all Me Z 2 and n 6 Z are summed 
over. The integrations can be extended to 9 € R if the 
summation over the M is restricted to a summation over 
a subset M. The selection of the subset fixes a gauge for 
M, which is the vector potential of the vortex density, 
iV = V x M. The summation over the subset of M then 
corresponds to a summation over the vortex density. 

The integration over the angles is most conveniently 
performed in Fourier space and results in the Hamilto- 
nian ([l2|). In the following we examine closer the vortex- 
charge interaction that reads on at an intermediate level 

H vc = 2mp J ^™-kik • M k . (A2) 
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Note that it is not possible to choose a Coulomb gauge 
with k-M k = due to the integer nature of M r . The sta- 
tistical weights exp(— Ti.) are invariant under a change of 
gauge M — > M+x due to the integer nature of the charge 
field n r and the gauge field \r- To transform expression 
( [A2] ) back into real space, let us look for simplicity at the 
ground state of Hamiltonian (Al) in the absence of fields 
and disorder, 



g(0) 
7 k 



2m. 



Ml 



(A3) 



But this is just the phase field generated by vortices, 
which can be written in real space (neglecting discretiza- 
tion effects) 



0<°> 



R 



R)JVi 



(A4) 



where uj(r — R) measures the angle enclosed between the 
vector (r — R) and an arbitrary reference di rect ion in the 
xy plane. In direct analogy the coupling ( |A2| ) becomes 
( |12cj ) in real space. 9^ and H vc do not depend on the 
choice of this reference direction for vortex configurations 
satisfying the neutrality (pf). Only these configurations 
have finite energy and contribute to the statistics. 



Then it is straightforward to perform the integration over 
the annulus 1 < |Ri — R2I < 1 + dl. After that step one 
can translate the vortex pair over the whole space ex- 
cluding the spheres occupied by the other vector pairs 
and obtains a screening of charges by vortices very simi- 
lar to the screening of vortices by vortices. Cj 



APPENDIX C: DICTIONARY 



For a convenient comparison of the present work with 
CO we include the following substitution list for symbols 
(restricted to the replica limit n — ► 0). 



CO symbol 
K 
K 
hp 
ft 

Vp 

V 
,,2 



symbol of present work : 
K, 

K — a T K 2 , 
n/%/2, 



(CI) 



V^pycUi 

47T 3 ^. 



Voth 



APPENDIX B: SCREENING OF CHARGES BY 
VORTICES 



Here we show how the angular interaction between vor- 
tices and charges induces a screening ( p^) of the interac- 
tion among charges when vortices are integrated out from 
the partition sum. Following KosterlitzEj and YoungEJ 
we integrate out every pair of vector vortices with radius 
less than 1 + dl. The contribution of a vortex pair with 
Ni + N2 = to the Hamiltonian TC VC is 

H vc = ipJ2 N ^ ' n r [w(Ri -r) -cj(R 2 - r)]. (Bl) 



Since we assume a diluted Coulomb gas, so that typically 
|Ri — r| 3> |Ri — R2I, one can expand the parts of the 
partition sum including this pair: 



exp (-« vc ) « 1 - — ^2 nr ' nr ' 

r,r' 

x [cj(Rx - r) - uj{R 2 - r)][w(R 1 - r') - uj(R 2 - r% (B2) 



In contrast to the screening effects within the vortices, 
the first nontrivial term in this expansion, which is even 
in the charge density and produces nonvanishing screen- 
ing, is of second order in the Hamiltonian. In the dilute 
limit we can expand 
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FIG. 1. Generic phase diagrams (schematically) for the 
special cases of (a) uniform fields for p > 4 and (b) random 
fields for p > 2^2- 




FIG. 2. This diagram indicates the location of the lines 
on the low-temperature side of which various char ges classes 
become relevant according to the scaling (24b). In the 
shaded area all charges are irrelevant. In the whole param- 
eter plane the most relevant charge class is either f (on the 
high-temperature side of the dashed line at K~ x = a ) or \\ 
(on the low-temperature side of the dashed line). 




FIG. 3. Nonreentrant phase diagram in the absence of 
fields according to Refs. |l9|-p3| and as it is reproduced by the 
flow equations fcj). 
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FIG. 4. Phase diagram in the absence of vortices. The 
transition lines between the phases LRO, QLRO, and QSRO 
depend only weakly on small fields Ho <§C 1. The transition 
between LRO and QSRO is almost horizontal for weak fields. 
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FIG. 5. Phase diagram for the model with a = a —a 
and uniform fields. It is obtained from a numerical integra- 
tion of the flow equations ( |49[ ) for two different values of Ho 
and 7 = 1.6. 
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